Outflows from dynamo-active protostellar accretion discs 
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Abstract. An axisymmetric model of a cool, dynamo-active accretion disc is ap- 
plied to protostellar discs. Thermally and magnetically driven outflows develop that 
are not collimated within 0.1 AU. In the presence of a central magnetic field from 
the protostar, accretion onto the protostar is highly episodic, which is in agreement 
with earlier work. 
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1. Introduction 

Outflows from protostellar discs have been observed in many star- 
forming regions. These outflows are generally believed to be driven 
by magneto-centrifugal acceleration from the accretion disc around the 
protostar. Yet, most theoretical approaches to understanding stellar 
outflows have ignored or at least greatly oversimplified both accretion 
disc physics and magnetic field evolution. Instead, most effort so far 
has focused around the physics of magnetised outflows and the colli- 
mation process assuming a fixed magnetic field. One of the perhaps 
most convincing simulations of the collimation of outflows have come 
from Ouyed et al. (1997) and Ouyed &: Pudritz (1997) who considered 
a disc corona with an imposed field aligned with the rotation axis of 
the disc and modelled the disc physics as a boundary condition. 

In recent years several groups have begun to approach the problem 
from different directions. One approach is to start from accretion disc 
physics and to extend the simulation domain in the vertical direction, 
so as to capture parts of the outflow acceleration region. The global 
simulations of Hawley (2000) are among the most impressive attempts 
to treat the problem as a whole. Here the magnetic field is generated 
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by the disc dynamo driven by the magneto-rotational instability and 
some form of outflow is indeed seen in the vicinity of the rotation 
axis. However, the dynamo still seems to be unable to produce a large 
scale poloidal field (which seems to be required for driving systematic 
outflows), and thermodynamics is not included to allow the disc to cool 
sufficiently to become geometrically thin. Other groups have focused 
on the dynamo aspect and have used thin disc theory as input to their 
models (Reyes-Ruiz Sz Stepinski 1995) or have studied the effect of 
vertical outflows on the disc dynamo (Bardou et al. 2001). 

In the present paper we review the approach of von Rekowski et al. 
(2003) which is conceptually closest to the dynamo approach mentioned 
above. The effect of the Lorentz force on the gas motions is taken into 
account. For a thin accretion disc, it is required that near-Keplerian 
shear velocity be the result of a balance mainly between centrifugal 
force and gravity. This is only possible if the disc is cool, and remains 
cool notwithstanding the adiabatic heating associated with the concen- 
tration of gas near the midplane by the vertical component of gravity 
during the disc formation. 

In order to stay conceptually as close as possible to the polytropic 
setup of Ouyed & Pudritz (1997) we adopt a bi-polytropic model with 
low entropy in the disc and high entropy in the exterior which we refer 
to as halo or corona (see the next section). 

Since we restrict ourselves to two-dimensional geometry, we model 
the disc dynamo using mean-field (a£l dynamo) theory. The appli- 
cability of this theory to protostellar accretion discs, where the tur- 
bulence (and the dynamo) is driven by the magneto-rotational (or 
Balbus-Hawley) instability, will also be discussed in the next section. 



2. Modelling a cool dynamo-active accretion disc 

Like Ouyed & Pudritz (1997) we start with an equilibrium corona, 
where we assume constant entropy and hydrostatic equilibrium accord- 
ing to CpT(r) = GM*/r. (Here, T is temperature, G is the gravitational 
constant, M* is the mass of the protostar, r is the spherical radius, and 
c p is the specific heat at constant pressure.) In order to make the disc 
cooler, we prescribe a geometrical region for the disc (see Fig. 1) and 
choose an entropy contrast between disc and corona such that the initial 
disc temperature is about 3 x 10 3 K in the bulk of the disc. Given the 
degree of simplification of our model, we have refrained from modelling 
the flaring of the disc. The low disc temperature corresponds to a high 
disc density of about 10~ 10 . . . 10~ 9 gem -3 . 
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Figure 1. General structure of the outflows typically obtained in our model. The 
cool, dense disc emits (a) a thermally driven wind (slow, hot, dense, magnetised, 
rotating), (b) a magneto-centrifugally driven outflow in a conical shell (faster, cooler, 
less dense, magnetised, quickly rotating), and (c) a thermally driven outflow near 
the axis (slow, hot, dense, weakly magnetised and weakly rotating). 

Hydrostatic equilibrium requires (tt- V)tt = — (l/p)Vp — V3>, where 
u = (0, u^,0) is the velocity field, p is gas pressure, p is gas density, 
and $ is the gravitational potential. Since specific entropy s is given, it 
is convenient to express (l/p)'Vp in terms of s and specific enthalpy h, 
using the relation V/i = (l//?)Vp + TVs (first law of thermodynamics) 
and h = c p T. In this way we obtain h by integrating the vertical balance 
equation, 

-d(h + ^)/dz + hd(s/c p )/dz = 0, (1) 

from large z (where h + $ = 0) down to z = 0. The initial rotation 
velocity, u v q, follows from the radial balance equation, 

-u^o/w = -d(h + <S>)/dw + hd{s/cp)/dw. (2) 

In the disc, h = c p T is small, so u^q is close to the Keplerian velocity, 
while the corona does not rotate initially, and is supported by the 
pressure gradient. Within the disc, we adjust density with a forcing 
term, so that a similar disc structure is maintained on reasonable time 
scales. 

For the disc dynamo, the most important parameter is the dynamo 
coefficient a in the mean-field induction equation. The a effect is an- 
tisymmetric about the midplane, restricted to the disc, and roughly a 
negative multiple of the kinetic helicity, so a is normally positive in 
the upper disc half. However, there is evidence that in accretion discs, 
where turbulence is driven by the magneto-rotational instability, a is 
negative in the upper disc half (Brandenburg et al. 1995, Ziegler &: 
Riidiger 2000). A simple theoretical explanation is given in Branden- 
burg (1998), and a similar result has been obtained by Riidiger k, Pipin 
(2000) who considered magnetically driven turbulence; see also Riidiger 
et al. (2001). 
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Is there any evidence that the aQ, dynamo picture is correct and 
applicable to accretion discs? One piece of evidence is that the specific 
dependence of the large scale field geometry on the boundary conditions 
is the same both in simulations of shearing sheet turbulence driven by 
the magneto-rotational instability and in aQ dynamo models in the 
same geometry. Since the dynamo-generated magnetic field is of large 
scale and since it extends to the boundaries, one must expect a de- 
tailed dependence of its geometry on the boundary conditions. In discs 
with negative a in the upper half, with (pseudo-)vacuum boundary 
conditions one gets an oscillatory large scale field whose horizontal com- 
ponents are symmetric about the midplane. With perfect-conductor or 
finite-conductivity boundary conditions, one gets a non-oscillatory field 
whose horizontal components are antisymmetric about the midplane. 
The same is seen both in three-dimensional simulations and in solutions 
of the a£l dynamo problem; see Brandenburg (1998). 

Bardou et al. (2001) found that vertical outflows do not change the 
symmetry or temporal behaviour of the dynamo-generated magnetic 
field. However, especially in the linear case for negative a in the upper 
half and finite-conductivity boundary conditions, a strong outflow tends 
to switch the dynamo off unless the dynamo number is large enough. 



3. Structured outflow and acceleration mechanism 

Figure 2, left panel, illustrates that an outflow develops that has a 
well-pronounced structure. Within a conical shell originating from the 
inner edge of the disc, the terminal outflow speed exceeds 500 km s" 1 , 
and temperature and density are lower than elsewhere. The inner cone 
around the axis is the hottest and densest region where the stellar wind 
speed reaches about 150 km s -1 . The wind that develops from the outer 
parts of the disc has intermediate values of the speed (see the sketch 
in Fig. 1). 

The structured outflow is driven by a combination of different pro- 
cesses. A significant amount of angular momentum is transported out- 
wards from the disc into the wind along magnetic field lines, especially 
along the strong lines within the conical shell (see Fig. 2, right panel). 
The magnetic field geometry is such that the angle between the rotation 
axis and the field lines threading the disc exceeds 30° at the disc surface, 
which is favourable for magneto-centrifugal acceleration (Blandford & 
Payne 1982). However, the Alfven surface is so close to the disc surface 
at the outer parts of the disc that acceleration there is mainly due to 
the pressure force. In the conical shell, however, the outflow is highly 
supersonic but sub-Alfvenic, with the Alfven radius a few times larger 
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Figure 2. Left: poloidal velocity vectors and poloidal magnetic field lines (white) 
superimposed on a colour scale representation of log 10 h. Specific enthalpy h is 
directly proportional to temperature T, and log 10 h — (—2, —1, 0, 1) corresponds to 
T w (3xl0 3 ,3xl0 4 ,3xl0 5 ,3xl0 6 )K. The black dashed line shows the fast magne- 
tosonic surface. The disc boundary is shown with a thin black line, the stellar surface 
is marked in red. The dynamo a coefficient is negative in the upper disc half, resulting 
in roughly dipolar magnetic symmetry. Averaged over times t « 897 . . . 906 days. 
Right: colour scale representation of the specific angular momentum, normalised 
by the maximum angular momentum in the disc, with poloidal magnetic field lines 
superimposed (white). The black solid line shows the Alfven surface, the white 
dashed line the sonic surface. Same model as left and averaged over same times. 



than the radius at the footpoint of the field lines at the disc surface. The 
lever arm of about 3 is sufficient for magneto-centrifugal acceleration 
to dominate (cf. Krasnopolsky et al. 1999). The inner cone has very 
low angular momentum and the magnetic lines are inclined to the 
axis by less than 30° so that the stellar wind can only be pressure 
driven. This picture is confirmed by looking at the ratio of magneto- 
centrifugal to pressure forces, which is much larger than unity in the 
conical shell. This is also the region where magnetic pressure produced 
by the toroidal field exceeds gas pressure, indicating that the outflow 
in the conical shell is confined by the toroidal field. 



4. Magnetic coupling between star and disc 



Observations of T Tauri stars suggest that the stellar spin-down can be 
connected with a circumstellar accretion disc. Magnetic fields are be- 
lieved to play an important role in the interaction between the star and 
its surrounding disc, resulting in particular in a spin-down of the star 
due to magnetic braking. Indeed, magnetic fields as strong as 10 3 Gauss 
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Figure 3. Colour scale representation of the density (bright colours indicate high 
values; dark colours, low values) with poloidal magnetic field lines superimposed 
(white) and the azimuthally integrated mass flux density shown with arrows (not 
shown in the disc). Left: at a time t w 150 days, when the star is magnetically 
connected to the disc. Right: at a time t ~ 156 days, when the star is disconnected 
from the disc. 



have been detected on T Tauri stars, and there is evidence for hot and 
cool spots on the stellar surface. Although the structure of the stellar 
magnetosphere is yet unknown, analytical models have been developed 
by Hartmann et al. (1994) and Shu et al. (1994) that assume that a 
stellar dipole threads a surrounding disc. In these static magnetospheric 
accretion models, disc matter is channelled along magnetospheric field 
lines to form a funnel-shaped polar accretion flow. 

More realistic time dependent numerical simulations by Hirose et al. 
(1997) and Goodson et al. (1997) show episodic reconnection behaviour 
between a stellar dipolar magnetosphere and an imposed disc field, 
resulting in episodic mass transfer from the disc to the star. These 
authors find inner stellar jets and outer disc winds that are both driven 
by magneto-centrifugal processes. 

Figure 3 shows our model with disc dynamo, where we add a stel- 
lar dipole that is anchored within the star, but can evolve with time 
outside the star (von Rekowski & Brandenburg 2003). For the strong 
surface stellar field of about 6kG considered here, the time dependent 
accretion is correlated with the magnetic star-disc coupling. When the 
star is connected to the disc, accretion is along field lines (see Fig. 3, 
left panel). When the star is disconnected from the disc, an enhanced 
wind carries away excess angular momentum (see Fig. 3, right panel). 
The sum of magnetic and accretion torques in the outward direction is 
positive at radii between star and disc, which means a stellar spin-down. 
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Figure 4- Colour scale representation of the normalised azimuthal velocity (top) 
and temperature (bottom) with poloidal magnetic field lines (white) and poloidal 
velocity vectors superimposed. Left: simulations of outflows from protostellar accre- 
tion discs by von Rekowski & Brandenburg (2003). Same model as in Fig. 3, but 
at t ~ 158 days. Collimation is due to boundary conditions. Right: simulations of 
outflows from accretion discs around black holes by Koide et al. (1998). 

5. Concluding remarks 

Our simulations show that a system consisting of a magnetically cou- 
pled star and dynamo-active accretion disc can maintain a clearly 
structured outflow. The outflow consists of a (hotter and denser) slow 
stellar wind and a faster outer disc wind that are both pressure driven, 
and a (cooler and less dense) fast magneto-centrifugally accelerated 
inner disc wind within a conical shell (see Fig. 4, left panels). In ad- 
dition, in transition periods, when the magnetic star-disc connectivity 
changes, there is a pressure driven (small azimuthal velocity and an- 
gular momentum, the Alfven surface close to the disc), hot and dense, 
but relatively fast outflow between the stellar and disc winds. 

Coexisting pressure and magneto-centrifugal driving mechanisms 
have also been found by Koide et al. (1998) who study the forma- 
tion of relativistic jets from accretion discs surrounding a black hole 
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by performing numerical simulations using a general relativistic mag- 
netohydro dynamic code. In their model, a global (initially vertical) 
magnetic field penetrates the accretion disc which then leads to non- 
steady accretion and ejection processes. According to their simulations, 
a jet is ejected from the close vicinity of the black hole. The jet has 
a two-layered shell structure consisting of a fast gas-pressure driven 
jet in the inner part and a slow magnetically driven jet in the outer 
part, both of which are collimated by the global poloidal magnetic field 
penetrating the disc (see Fig. 4, right panels). The fast jet is due to a 
strong pressure increase due to shock formation at the inner edge of 
the disc, caused by a fast accretion flow. 

Our relatively fast pressure driven flow might be comparable with 
the fast jet. However, contrary to the slow jet, our magneto-centrifugally 
accelerated inner disc wind is even faster. 
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